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Abstract 

The Ultrarelativistic Quantum Molecular Dynamics (UrQMD) model, a microscopic transport 
model, is used to study the directed and elliptic collective flows and the nuclear stopping in Au+Au 
collisions at incident energies covered by INDRA and lower-energy FOPI experiments. It is seen 
clearly that these observables are sensitive to both, the potential terms (including iso-scalar and 
iso- vector parts as well as the momentum dependent term) in the equation of state (EoS) and 
the collision term (including the Pauli-blo eking and the medium-modified nucleon-nucleon elastic 
cross section (NNECS)). The momentum modifications of both, the mean-field potentials and the 
density dependent NNECS, are found to be sensitive to the collectivity of heavy-ion collisions. At 
INDRA energies (< 150 MeV/nucleon), the dynamic transport with a soft EoS with momentum 
dependence and with the momentum-modified density-dependent NNECS describes the directed 
flow exhibited by hydrogen isotopes (Z = 1) emitted at mid-rapidity fairly well. 

PACS numbers: 25.70.-z,24.10.-i,25.75.Ld 
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I. INTRODUCTION 



Nuclear reactions at intermediate energies as, e.g., provided by the heavy- ion synchrotron 
SIS at GSI have been investigated for several decades but the physics governing their mecha- 
nisms are not yet thoroughly understood and further investigations seem necessary. Recently 
the two experimental collaborations INDRA and FOPI have published a series of systematic 
observations made in experiments at the GSI laboratory regarding several physical quantities 
such as the collective flows, the nuclear stopping, and the light cluster and produced-particle, 
i.e. mostly pion, production [l-4]. Comparisons with model calculations (mainly performed 
with the Quantum Molecular Dynamics (QMD) model), revealed some discrepancies which 
deserve deeper investigations, mainly concerning dynamical observables. The excitation 
functions of the directed and elliptic collective flows and the nuclear stopping cannot be 
satisfactorily described with the same equation of state (EoS) of nuclear matter (see e.g. 

Important progress in nuclear physics at intermediate energies was made recently by 
reaching consensus on the soft nature of the EoS of nuclear matter j^, Q-[lQ]. However, 
according to the theory of quantum hydrodynamics (QHD), both the mean field and the 
two-body collisions have the same origin, the effective Lagrangian density, based on which 
the medium modifications of both the collisions and the mean- field transport between the col- 



lisions should be considered self-consistently 
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3- 



More explicitly, the medium-modified 



terms should take the density dependence, the isospin asymmetry, and the momentum con- 
straints into account. In past investigations, these medium effects were considered mainly 
in the mean-field part but usually not in the collision part, leading to an obvious lack of 
self-consistency. Hence, conclusions based on these treatments are not fully reliable, and 
it is easy to understand why new difficulties arose in the recent comparisons of data with 
model calculations. 

In the lower range of SIS energies, the collision rate is known to increase with beam 
energy and the interplay between mean field and two-body collisions leads to a colorful 



phenomenology o 



collective flows 



' the nuclear dynamics as, e.g., apparent in the excitation functions of 



16 



3Hl8|: with the increase of beam energy beyond the INDRA range (< 
150 MeV/nucleon), the slope, near mid-rapidity and with respect to rapidity, of the directed 
flow of free protons or light clusters changes sign from negative to positive, while the value 



of the elliptic flow at mid-rapidity changes from positive (in-plane) to negative (squeeze- 
out). These transitions of the directed and elliptic flows imply clearly that the strength of 
two-body collisions and the nuclear stopping power increase, which has been observed in 



recent INDRA experiments as well 



19j. Theoretically, in addition to the well-known 



density- dependent modifications of binary cross sections, the momentum constraints will 
have to be considered seriously in order to describe these data systematically 2CH31]. 

In our previous work, we have attempted to build up a microscopic transport model in 
which the dynamic process of heavy ion collisions (HICs) is described in a more consistent 
and complete manner. The Ultrarelativistic Quantum Molecular Dynamics (UrQMD) model 



is adopted as a microscopic transport basis 32l434j . In Refs. 



35 



both the mean-field 



potentials and the relativistic effects on the relative distance and momentum in potentials 



are treated comprehensively. In Refs. 



26 



the medium modifications of the nucleon- 



nucleon elastic cross sections (NNECS) are further taken into account. With the so improved 
version of the UrQMD model, several experimental observables at SIS and AGS energies were 



successfully described 
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42j. With the further consideration of the "pre- formed" hadron 



potentials, the well-known HBT time-related puzzle throughout the beam energy range 
from SIS up to RHIC is well understood 43|, |44j]. It is also found, however, that further 
investigation is still needed at the lower SIS energies [40|. This will become evident in this 
work from the study of collective flows and nuclear stopping. 

The paper is arranged as follows: Section II presents the new updates in the UrQMD 
transport model. Results for collective flows and nuclear stopping in the energy range covered 
by INDRA and FOPI experiments are then shown in Section III. Finally, a conclusion and 
an outlook are given in Section IV. 



II. URQMD TRANSPORT MODEL AND UPDATES 



The UrQMD model is based on analogous principles as the QMD model 



46| and the 



Relativistic Quantum Molecular Dynamics (RQMD) model j47J. The first forma 



32 



version 



33. 



48]. 



(verl.O) of the UrQMD transport model was published at the end of last century 
Since then, a large number of successful theoretical analyses, predictions, and comparisons 
with data have been accomplished with this transport model for pp, pA and AA reactions 



and over a large range of beam energies, essentially from SIS over AGS, SPS, RHIC, up to 



4 



LHC energies 



49] 



A. Potential updates 

In the UrQMD, similar to the QMD, hadrons are represented by Gaussian wave packets 
in phase space. After the initialization of projectile and target nuclei for which the Hard- 
Sphere (H-S) or the Woods-Saxon (W-S) mode may be selected, the phase space of hadron 
i is propagated according to Hamilton's equations of motion: 

. dH . dH 

r 4 = 5—, and p; = -— . (1) 

opt dvi 

Here r and p are the coordinate and momentum of hadron i. The Hamiltonian H consists 

of the kinetic energy T and the effective two-body interaction potential energy V, 

H = T + V, (2) 

with 



T = - m) = J2W m " + Pi - ( 3 ) 



In the standard version of the UrQMD model [32j, |33j, |48|], the potential energies include the 



two-body and three-body Skyrme-, and the Yukawa-, Coulomb-, and Pauli-terms as a base, 

V = Kg + Vg + V Yuk + V Cou + V Pau . (4) 

The three-body terms of the Skyrme force can be approximately written in the form of two- 
body interactions, the so-called density dependent terms. Only the Coulomb force between 
charged baryons is considered in the default version. However, when also considered for 



charged mesons it was found indispensable for reproducing certain observables 35l. l38|. 

The single particle potential follows from U = 5V/5f, where / is the phase-space Wigner 
distribution function which reads as 

/(r, p) = Mr, P) = E 7i^ (r ^ l)2/2L2e " (P ' P!)2 - 2L2/?i2 ' ( 5 ) 

i i ^ ' 

where L is the width parameter of the wave packet. It is known that the width affects 

n 

the stability of initialized nuclei and the production of clusters at later reaction stages 50J. 



Empirically, the value of the wave packet width should grow from about 1 fm to 2 fm as the 
total mass of the colliding system increases from light to heavy. As an example, for Au+Au 
collisions a value of L = 2 fm is chosen while for Ca+Ca and lighter nuclei L = 1 fm is found 
more appropriate. 

Recently, in order to permit more successful applications in the intermediate energy 
regime (0.1 A < E lab < 2 A GeV), more potential terms have been incorporated into the 
UrQMD 37|. The first addition is the density-dependent average symmetry potential energy 
term V^S 2 , where 5 = (p n — p p )/po is the isospin-asymmetry. In this work, the form 



ypot 
sym 



(S 



e F , 



(6) 



introduced in Ref. 



51| is used. Here So is the symmetry energy at the normal nuclear density 



Po, ep is the Fermi kinetic energy at normal nuclear density, u = p/po is the reduced nuclear 
density, and 7 is the strength of the density dependence of the symmetry potential. In this 
work, the value Sq = 34 MeV is chosen. The symmetry potential is important for isospin- 
asymmetric reactions at intermediate and low energies since the uncertainties existing in 
the symmetry energy based on various theories were found to be large at densities away 



from the normal density 



52] . Here, a soft (7 = 0.5) density dependent symmetry potential 



is adopted. Recent investigations indicate that the density dependence of the symmetry 
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energy at subnormal densities might be soft [29 
confirmed by further studies. 

The second addition is the momentum-dependent term 



54J. However, this will still have to be 



55] 



Umd = tmd In 2 [1 + d m d(Pi ~ Pjf]pi/ ' (h 



0; 



(7) 



where t m d = 1.57 MeV and a m d = 500 c 2 /GeV 2 are selected as in previous calculations 
within the QMD family. It has to be noted, however, that this form may also have to be 
refined in order to obtain improved fits of the real part of the optical potential js), [ssj ] . 

At beam energies higher than the AGS regime, potentials for "pre-formed" particles 
(string fragments) and relativistic effects regarding the relative distance and the relative 
momentum em ploy ed in the two-body potentials (Lorentz transformation) are, furthermore, 



considered 
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43]. They can be neglected in this work here. 



B. New Pauli-blocking treatment for two-body collision 



Regarding the collision term, it is well-known that the rate of collisions is reduced by the 
quantum Pauli-blocking. However, in the default version of the UrQMD, for simplicity, the 
Pauli-blocking criteria are realized with an algorithm based on fitted parameters: for the 
same type of particles, the binary scattering will be blocked under the form 



V < 4>, (8) 

where 

0- r ,i) 2 (P-Pi) 2 -£ 2 

ip = e 4L ' 2 e ^ , (9) 

i 

and 



V = a fit + b f itQ (10) 

with 

e (11) 

i 

The parameters afu and bfu in Eq. f JTOj) are set to be 1.49641 and 0.208736, respectively. It 
will be shown below that, for descriptions of experimental observables at the lower energies, 
the default Pauli-blocking treatment is not precise enough. Therefore, an updated new 
version of the Pauli-blocking routine is introduced here. For the same type of particles and 
for each collision, firstly, the phase space densities in the final states are determined in order 
to assure that they are in agreement with the Pauli principle and, secondly, the following 
two criteria are considered at the same time: 

—r% ■ — V % > (2s + 1) • (-) 3 (12) 

and 

Puock = i-a -mi- fj)<t (13) 

On the left-hand side of Eq. (1121) . the and Pij are the relative distance and momentum of 
the two particles at the final states i and j. The factor (2s + 1) on the right side of Eq. ( TT2"j) 
denotes the summation of the spins of the two particles. In Eq. ( 1131) . £ represents a random 
number between and 1. If one of the criteria is not fulfilled the collision is not allowed and 
the two particles remain with their original momenta. 



Set 


A 


c 


FU1 


1/3 


0.54568 


FU2 


1/4 


0.54568 


FU3 


1/6 


1/3 



TABLE I. The three parameter sets FU1, FU2, and FU3 used for the density-dependent correction 



factor F u of NNECS. 



C. Medium modifications and momentum dependence of NNECS 



py the nuclear medium, according 
. As in previous work, the in- 



It is also known that the cross sections will be modifiec 
to the QHD theory (see, e.g., Refs. \m\m, hi M M |57 
medium NNECS a* d are treated to be factorized as the product of a medium correction 
factor F(u,a,p) and the free NNECS crf/" ee 26]. For the inelastic channels, we still use 
the experimental free-space cross sections cx{[ ee . It is believed that this assumption does 
not have strong influence on our present study at low SIS energies. Therefore, the total 
two-body scattering cross section of nucleons, a^ ot , will be modified according to 



(14) 



As for the medium correction factor F(u,a,p), it is proportional to both the isospin-scalar 
density effect F u and the isospin- vector mass-splitting effect F a . Studies of the isospin- 
related split ting effect on NNECS which is represented by the F Q factor have been reported 



in Refs. 



26 



Furthermore, the factors F u and F a should be functions of the relative momentum pnn of 



the two colliding particles in the NN center-of-mass system. In Ref. 
as 



26] , they are formulated 



/o 



-fo 



l+(pjviv/po) 

The factor F u can be expressed as 



- + fo 



Pnn > lGeV/c 
Pnn < lGeV/c. 



(15) 



F u = A + (l-A)exp[-u/C]- 



(16) 



S 



Set 


fo 


Po [GeV c- 1 ] 


K 


FP1 


1 


0.425 


5 


FP2 


1 


0.225 


3 


FP3 


1 


0.625 


8 


no pnn limit 


F(u) 


/ 


/ 



TABLE II. The three parameter sets FP1, FP2, and FP3 used for describing the momentum 
dependence of F u . The fourth case, without a pnn limit, is obtained by setting /o equal to F(u) 
in Eq. (fT5l) . 



Here ( and A are parameters which determine the density dependence of the cross sections. 
In this work, we select several parameter sets which are shown in Table UJ and illustrated 
in the left panel of Fig. [TJ The reduction of the NNECS as a function of density becomes 
increasingly more pronounced as the parameterization is changed from FU1 to FU3. At 
the reduced density u = 2, e.g., the values of FU1, FU2, and FU3 are 0.35, 0.27, and 0.17, 
respectively. We note here that the density dependence of the FU1 parameterization is in 



cjualitative agreement with previous work based on the Dirac-Brueckner approach 30 



31 



59]. However, in our previous investigations of the NNECS, based on the effective Lagrangian 



of density-dependent relativistic hadron theory in which the a, u, p and 5 [ao(980)] mesons 
are included [23], it was shown that especially the neutron-proton cross sections <J* l np might 
be largely reduced in the neutron-rich nuclear medium; the corresponding reduction factor 
might be as low as ~ 0.1 at u = 2. Therefore, the other parameter sets FU2 and FU3 (Table 
UJ are still to be considered reasonable assumptions. 

The parameters /o, po and k in Eq. ( jT5|) can be varied in order to obtain various mo- 
mentum dependences of F u . In this work, we select several parameter sets which are shown 
in Table [Til The corresponding F£ functions are illustrated in the right panel of Fig. [1] for 
FU1 at a reduced density u = 2. The FP1 set was used in our previous work in which the 



26]. The parameterizations FP2 or 



medium modifications of cross sections were considered 
FP3 result in a rise of the correction factor F at smaller or at larger paw than obtained 
with FP1. It illustrates the present uncertainty associated with treating the momentum 
dependence of the density-dependent cross sections. With a certain set of isospin dependent 
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FIG. 1. (Color online) Correction factor F obtained with the parameterizations FU1, FU2, and 
FU3 given in Table U (left panel) and the momentum dependence obtained with the four options 
FP1, FP2, FP3, and "no pnn limit" given in Table HI1 for the example of FU1 at u = 2 (right 
panel). 



EoS, the NNECS might even be enhanced at large momenta, in comparison to the cross 
sections at free space. It arises from the differences between the isoscalar and isovector 



channels 



151 ]. An indication of this enhancement was obtained from a recent calculation of 



flow observables by Zhang et al. 



20j. But it will not be discussed in the current work since 



the effect is related to collisions at large Pat at and it becomes increasingly important for 
heavy-ion collisions at higher beam energies. And, for completeness, the case without any 
momentum constraint on F u ("no Pnn limit") is listed in Table HT1 
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III. CALCULATIONS AND OBSERVABLES 



A. Colliding system and the after-burner 

About 120 thousand events of 197 Au+ 197 Au collisions for each of the energies Ei a t> = 40, 
50, 60, 80, 100, 150, and 400 MeV/nucleon are calculated randomly within the impact 
parameter region 0-7.5 fm. For each beam energy, the calculations are divided into three 
groups according to impact parameter representing the central (0-2 fm), the semi-central (2- 
5.5 fm), and the semi-peripheral (5.5-7.5 fm) group of collisions. Fourteen parameterizations 
of the UrQMD transport model differing in the treatments of the Pauli-blocking of two- 
body collisions (P-B), of the potential terms (EoS), of the medium- modified NNECS with 
the density-dependent factor F u and the momentum modification F p of F u , and of the 
initialization (ini.) are shown in Table HTT1 and will be discussed in the following section. The 
UrQMD transport program stops and produces the output after a collision time of 150 fm/c. 
Several more outputs are produced, as the calculations evolve, at different reaction times 
before the stop time. For each output, a conventional phase-space coalescence model \§0\ 
with two parameters is used to construct clusters. Nucleons with relative momenta smaller 
than Po and relative distances smaller than Rq are considered to belong to one cluster. In 
this work, P and Rq are chosen to be 0.2 GeV/c and 2.8 fm, respectively. It is evident 
and well-known that the criteria chosen for constructing clusters affect their yields and, 



consequently, also related observables [61 ]. 



B. Collective flows and the stopping vartl 

In this work, we will focus on the parameters describing the strength and orientation of 
directed and elliptic collective flows defined as 

Vl =< cos((/> - <S> RP ) >=< — >, (17) 

Pt 

and 

2 2 

v 2 =< cos[2(0 - $ flP )] >=< ^—fit > . (i 8 ) 

Pt 

Here denotes the azimuthal angle of the considered outgoing particle; ^/rp is the azimuthal 
angle of the reaction plane which is pre-configured in the calculations, ^rp = 0, but has to be 
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UrQMD-II 
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Cascade 


FU1 


FP1 


W-S 
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UrQMD-VI 
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SM 


no Collisions 


no Collisions 


H-S 


UrQMD-VII 


new 


SM 


FU1 


FP3 


H-S 


UrQMD-VI I I 


new 


SM 


FU1 


no pnn limit 


H-S 


UrQMD-IX 


new 


S 


FU1 


FP1 


H-S 


UrQMD-X 


new 


S 


FU1 


no Pnn limit 


H-S 


UrQMD-XI 


new 


S 


FU2 


no pnn limit 


H-S 


UrQMD-XII 


new 


S 


FU3 


no pnn limit 


H-S 


UrQMD-XI I I 


new 


SM 


FU3 


FP1 
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TABLE III. Fourteen parameterizations of the UrQMD transport model differing in the treatments 
of the Pauli-blo eking of two-body collisions (P-B), of the potential terms (EoS), of the medium- 
modified NNECS with the density-dependent factor F u and the momentum modification F p of F u , 
and of the initialization (ini.). See text for details. 

determined in experiments; p x and p y are the two components of the transverse momentum 
Pt = \fv1 + V 2 y - The angular brackets denote an average over all considered particles from 
all events. 

The nuclear stopping is described with the quantity vartl [62] denned as 

vart l = (19) 

LdN/dy z 
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Here 

Tfw 2 N ) 

?dN/dy x , z =< y x , z >= v Ar , (20) 

where T dN / dyx and V dN / dyz are the variances of the rapidity distributions of fragments in 
the x and 2 directions, respectively. They are obtained as weighted averages of y 2 over the 
considered ranges in rapidity with N Vx and N yz denoting the yields of fragments in each of 
the y x and y z rapidity bins. 

Directed and elliptic flows and the nuclear stopping quantity vartl were measured by 
the INDRA and FOPI collaborations at GSI/Darmstadt In the following, the pre- 

dictions of the UrQMD model for these observables will be analyzed, while the directed 
flow is chosen for the comparison with the experimental data at INDRA energies (Ei a b < 
150 MeV/nucleon). 



IV. RESULTS 

A. Pauli-Blocking effects on flows 

As stated before, the Pauli-blocking treated in the default UrQMD model is an effective 
one. The advantage of this treatment is the high-speed to calculate it. However, it is found 
that it is not quite suitable for nuclear reactions at low SIS energies. In Fig. [2] we present 
the reduced longitudinal rapidity distribution of directed (t>i, top plots) and elliptic (v 2 , 
bottom plots) flows of free protons from semi-peripheral (b=5.5-7.5 fm) Au+Au collisions 
at Eiab = 40, 80, and 400 MeV/nucleon. The calculations with the default treatment of 
Pauli-blocking used in UrQMD-0 and with the new treatment incorporated in UrQMD-II 
(the only difference between the two versions, cf. Table IIIII) are shown in Fig. [2j 

Considerable differences exist in the rapidity distribution of both flow parameters at the 
lower SIS energies E^b = 40 and 80 MeV/nucleon. With the default treatment, the slope of 
v\ is significantly positive at mid-rapidity and the value of v-i is rather small. Experimentally 
it is known that the slope of the directed flow of Z = 1 particles from semi-central/peripheral 
collisions of heavy systems at 40 MeV/nucleon is slightly negative at mid-rapidity and that 
the elliptic flow is positive with values near v<i = 0.08 [1], values much closer to the results 
obtained with the new treatment of Pauli-blocking (Fig. |2j left panels). It is, therefore, 
concluded that this modification of the UrQMD is more suitable for this energy range, 
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FIG. 2. (Color online) Parameters v\ of directed flow (top panels) and V2 of elliptic flow (bot- 
tom panels) for free protons from semi-peripheral (b = 5.5 7.5 fm) 197 Au+ 197 Au collisions at 

E lah = 40 (left), 80 (middle), and 400 MeV/nucleon (right), calculated with the versions UrQMD-0 
(squares) and UrQMD-II (circles), as a function of the reduced rapidity y z /y pTO - 

and the new treatment of Pauli-blocking is chosen for the following investigations. The 
differences start to disappear at the higher energy 400 MeV/nucleon (right panels of Fig. 

B. Causes of collective phenomena 

In order to understand how the initialization, the mean-field potentials and collisions 
influence collective flows, calculations were performed with the versions UrQMD-II, UrQMD- 
III, UrQMD-IV, UrQMD-V, and UrQMD-VI (Fig. [3}. Several interesting phenomena are 
observed: comparing first the results obtained with UrQMD-II and UrQMD-III which differ 
only in the initialization, small differences are seen but the directed flow is slightly less 
positive at 40 MeV/nucleon and slightly more positive at 80 MeV/nucleon with the Hard- 
Sphere treatment of UrQMD-II than with the Wood-Saxon initialization of UrQMD-III. 
This is due to the fact that at 40 (80) MeV/nucleon the net-contribution of both mean- 
field potentials and two-body collisions is more attractive (positive) and the Hard-Sphere 
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FIG. 3. (Color online) Parameters v\ of directed flow (top panels) and V2 of elliptic flow (bottom 
panels) for free protons from semi-peripheral (b = 5.5 — 7.5 fm) 197 Au+ 197 Au collisions at = 40 
(left) and 80 MeV/nucleon (right) as a function of the reduced rapidity y z jy vm - The calculations 
performed with versions UrQMD-II, UrQMD-III, UrQMD-IV, UrQMD-V, and UrQMD-VI (see 
Table IIIip are distinguished by their symbols as indicated. 
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initialization may lead to stronger attraction (repulsion). Actually, this effect of different 



initializations was seen before in calculations with the IQMD model [50j . Since its effect on 
final flows is very small, mainly the Hard-Sphere initialization will be used in the following 
calculations. It is the default initialization in the UrQMD model. 

Secondly, if the mean-field potentials are ignored as in the cascade-type version UrQMD- 
IV, a large positive slope of v± is seen at both beam energies. It implies that the net 
contribution of the collision term to the directed flow is repulsive. The elliptic flow has 
almost vanished which should be due to both, a relatively weak contribution of the collision 
term and the almost isotropic differential NNECS. On the other hand, if the collision term 
is switched off and only the mean-field potentials are active as in UrQMD- VI, a strong 
negative directed flow is seen at both beam energies. The elliptic flow v% is largest in this 
case, indicating that the mean-field potentials contribute to the collective flows much more 
strongly than the collision term at lower SIS energies, a fact that is well-known. It is also 
found that the net contributions of the mean-field and of two-body collisions to the directed 
flow are attractive and repulsive, respectively. The positive elliptic flow induced by the 
mean-field is caused by the in-plane geometry at early reaction times. 

The isospin-scalar and isospin-vector related parts of the potentials may yield different 
contributions to the final flows. This was tested with the results of UrQMD-II, UrQMD-IV, 
and UrQMD-V in which a soft-EoS with momentum dependence (SM-EoS, the corresponding 
incompressibility K = 200 MeV) is considered, not considered, and partially considered, 
respectively, while the same collision term is adopted. It is found that the iso-scalar part of 
the potentials produces a strong attractive force (comparison of UrQMD-IV and UrQMD- 
V) while the iso-vector related part gives a relatively weak repulsive force to the directed 
flow (see from UrQMD-V to UrQMD-II). The repulsion caused by the iso-vector potentials 
(symmetry and Coulomb potentials) is obvious since the Coulomb potential of protons is 
always repulsive and much stronger than the symmetry potential, regardless of whether 
the total effect of the symmetry term is attractive or repulsive. The latter depends on the 
isospin-asymmetry of the nuclear medium during the whole dynamic process of HICs. 

Summarizing this analysis, one finds clearly that 1) the balance of contributions between 
mean-field potentials and the collision term determines the vanishing of the directed flow at 
low SIS energies with roughly equal strength; 2) in the mean-field potentials, the iso-scalar 
part behaves as an attractive effect which is due to the overall rather low densities of heavy- 

16 



ion collisions at such low beam energies. The iso-vector part is repulsive which is mainly 
due to the mutual Coulomb repulsion of the protons. 

Another set of comparisons of flow results, calculated with versions UrQMD-I, UrQMD-II, 
UrQMD-VIII, and UrQMD-IX is shown in Fig.HJ Two beam energies 40 and 400 MeV/nucleon 
are selected. The comparison of UrQMD-I, UrQMD-II, UrQMD-VIII shows the effects on 
the flow parameters resulting from different momentum-modified NNECS while the com- 
parison of UrQMD-II and UrQMD-IX shows those of the momentum-dependent terms in 
the potentials (with or without momentum dependence). It is, first of all, noticed that 
both effects, the momentum modification of the density-dependent NNECS in the collision 
term and the momentum dependence of the potentials affect both flow variables v\ and 
Secondly, when studied in more detail, these flow effects are found to behave differently at 
different beam energies: (1) regarding the contribution of medium- modified cross sections, 
we see that for the directed flow at 40 MeV/nucleon the difference is large between the 
UrQMD-I and UrQMD-II results, while at 400 MeV/nucleon the largest difference exists 
between the UrQMD-I and UrQMD-VIII results. It is obviously caused by the difference of 
the momentum modifications of cross sections (FP1, FP2, and no Pnn limit) as they appear 
at different energies; (2) for the elliptic flow at 40 MeV/nucleon, a difference is seen between 
the UrQMD-VIII and UrQMD-IX results, while this is no longer the case at 400 MeV. 
It implies that, at 40 MeV/nucleon, the momentum-dependent term in the potentials is 
more important than the momentum modification of the density dependent NNECS. They 
become equally important at 400 MeV/nucleon. 

A rough comparison of the calculations with INDRA and FOPI data for Z = 1 parti- 
cles from the same system shows that the effect of the momentum-dependent terms in the 
potentials is repulsive, i.e. not suitable for describing the measured negative slope of the 
V\ flow and the largely positive v% (in-plane flow) at mid-rapidity and at the beam energy 
40 MeV/nucleon. At 400 MeV/nucleon, the consideration of the momentum dependence 
in both, the mean-field and collision terms, is important to describe the measured largely 
positive slope of V\ and largely negative v 2 (squeeze-out) at mid-rapidity. Therefore, the 
uncertainties existing in the momentum modification of NNECS will largely influence the 
excitation functions of flows, even only at low SIS energies. 

The yield distributions of free protons as functions of the reduced longitudinal {y z /y W o) 
and transverse (y x /y P To) rapidities from central Au+Au collisions at = 80 MeV/nucleon 
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FIG. 4. (Color online) Parameters v\ of directed flow (top panels) and V2 of elliptic flow (bottom 
panels) for free protons from semi-peripheral (b = 5.5 — 7.5 fm) 197 Au+ 197 Au collisions at = 40 
(left) and 400 MeV/nucleon (right) as functions of the reduced rapidity y 2 /y P ro- The calculations 
performed with versions UrQMD-I, UrQMD-II, UrQMD-VIII, and UrQMD-IX (see Table [III]) are 
distinguished by their symbols as indicated. 
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FIG. 5. (Color online) Yield distributions of free protons as functions of the reduced longitu- 
dinal (y 2 /y P ro, squares) and transverse {y x /Vpro-, circles) rapidities from central (b = — 2 fm) 
197 Au+ 197 Au collisions at E Xsh = 80 MeV/nucleon. Calculations with UrQMD-0, UrQMD-I, 
UrQMD-II, and UrQMD-VIII (see Table HIT|) are shown in the top-left, bottom-left, top-right, and 
bottom-right panels, respectively. The corresponding values obtained for the stopping observable 
vartl for free protons are also indicated in each panel. 



are shown in Fig. Calculations with the default Pauli-blocking treatment in UrQMD-0 
(top-left) are compared to those with the new Pauli-blocking treatment combined with the 
same (UrQMD-II, top-right) or alternative momentum modifications of NNECS (UrQMD-I, 
bottom- left and UrQMD-VIII, bottom-right). The corresponding values of vartl, for the full 
range of rapidities, are also given for each of the UrQMD versions. We first notice that, in 
the UrQMD-0 mode, the relative maximum of the proton yield is higher at mid- longitudinal- 
rapidity than at mid-transverse-rapidity, which leads to vartl > 1. In contrast, calculations 
with the new treatment of the Pauli-blocking, despite of uncertainties in the cross sections, 
the values of vartl are always smaller than unity, in line with the experimental data. 

Secondly, with decreasing cross sections (the reduction of NNECS due to the momentum 
modifications increases successively going from UrQMD-I to UrQMD-II and UrQMD-VIII), 
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;he calculated values of vartl become smaller and approach the experimental value 0.80±0.03 
6] measured with INDRA at = 80 MeV/nucleon (note, however, that one can not 
compare directly with data since the experimentally selected clusters (Z < 3) and rapidity 
bins are different from the current calculations performed for Z = 1 particles). Together 
with the flow results shown in Fig. HI these comparisons clearly show that the momentum 
modification of NNECS are important for describing the collective flows and the nuclear 
stopping at the low SIS energies at the same time. Thirdly, in the analyses related to Figs. 
131 HI and [51 one finds that the results for flows and for the nuclear stopping always follow 
in the same order when different treatments of the mean-field and the collisions terms are 
chosen. It, therefore, seems sufficient to select one of the observables, the slope of the 
directed flow at mid-rapidity, for the following more precise comparison with experimental 
data. 



C. Comparison of the excitation function of the v% slope with INDRA data 

A modification of the density dependence might be considered as an additional alternative 
to the momentum modification of NNECS in order to achieve their stronger reduction in 
heavy-ion collisions at lower energies and a weaker reduction at higher energies. In the 
left half of Fig. El we show values of the slope of directed flow at mid-rapidity (\y z /y pIO \ < 
0.4) for Z — 1 particles from semi-central (b = 2 — 5.5 fm) Au+Au collisions at = 
40 MeV/nucleon. Calculations with UrQMD-VI, UrQMD-X, UrQMD-XI, and UrQMD-XII 
are compared to the INDRA data. When the collision term is turned off (UrQMD-VI), the 
absolute value of the v i slope is larger than the data, indicating that collisions start to become 
important for heavy-ion collisions at such low beam energies. When the collision term is 
present and larger reductions of NNECS are considered by choosing UrQMD-X, UrQMD- 
XI, and UrQMD-XII (with increasingly larger density dependence but without momentum 
modification), the obtained v± slopes clearly approach the measured INDRA value. 

On the right side of Fig. El again values of the slope of directed flow at mid-rapidity 
for Z = 1 particles from semi-central Au+Au collisions but for an incident energy Ei a b = 
150 MeV/nucleon are presented. Calculations with UrQMD-I, UrQMD-II, UrQMD-VII, and 
UrQMD-X are compared to the INDRA data. Obviously, with the increasingly weaker mo- 
mentum modification of NNECS (as realized with UrQMD-I, UrQMD-II, and UrQMD-VII 
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FIG. 6. (Color online) Values of the slope of directed flow at mid-rapidity (|y z /y P ro| < 0.4) for 
Z = 1 particles from semi-central (6 = 2—5.5 fm) 197 Au+ 197 Au collisions at = 40 MeV/nucleon 
(negative values, left) and 150 MeV/nucleon (positive values, right). The INDRA results at both 
beam energies are compared with UrQMD calculations performed with the indicated versions (see 
Table ITlTl) . all represented by the bars with different shadings. 



in that sequence), the large reduction induced by the density dependence of the NNECS is 
less compensated, leading to a weaker repulsion and, hence, to the decreasing value of the 
V\ slope shown in the figure. Alternatively, if we switch off both the momentum-dependent 
term in the potentials and the momentum modification of NNECS in the collision term, re- 
alized with the UrQMD-X calculation, the value of the v\ slope is pushed down further and 
approaches the INDRA data. As a result of these tests, one may reach a deeper understand- 
ing of the difficulties in theoretically describing experimental data and of the importance 
of treating self-consistently the dynamic process of heavy-ion collisions. Furthermore, it is 
evident from Fig. that even the version UrQMD-X is still not good enough to permit the 
description of the directed flow over the whole range of INDRA energies. Further refinements 
and tests will be needed. 

The time evolution of the slope of directed flow at mid-rapidity for Z — 1 particles from 
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FIG. 7. (Color online) Time evolution of the slope of directed flow at mid-rapidity (|y 2 /y P ro| < 
0.4) for Z = 1 particles from semi-central (6 = 2 — 5.5 fm) 197 Au+ 197 Au collisions at = 
40 MeV/nucleon (top) and 150 MeV/nucleon (bottom). Calculations with UrQMD-X and UrQMD- 
XIII (see Table IIIip are represented by the squares and circles, respectively. The INDRA results 
for both beam energies are indicated by the stars positioned at the UrQMD stop time of 150 fm/c. 
The dashed horizontal lines representing zero slopes are included for comparison. 

semi-central Au+Au collisions at = 40 (top) and 150 MeV/nucleon (bottom) are shown 
in Fig. [71 The versions chosen for the comparison are UrQMD-X, closest to the INDRA data 
at 150 MeV/nucleon (Fig. ED, and UrQMD-XIII. The latter differs from UrQMD-XII, closest 
to the data at 40 MeV/nucleon, in that the SM-EoS and the FP1 momentum modification 



of NNECS (the standard in previous investigations 26|) are adopted and from UrQMD-X 

22 




FIG. 8. (Color online) Excitation function of the slope of directed flow at mid-rapidity (|y z /y P ro| < 
0.4) for Z = 1 particles from semi-central (b = 2 — 5.5 fm) 197 Au+ 197 Au collisions in the range 
of INDRA energies from 40 to 150 MeV/nucleon. Calculations with UrQMD-X (dashed) and 
UrQMD-XIII (full line) are compared with the INDRA data (stars). The dotted horizontal line 
representing the zero slope is included for comparison. 

by the stronger density dependence FU3 for NNECS (same as in UrQMD-XII). At about 
t < 30 fm/c, the pre-equilibrated protons are emitted with a small positive flow due to 
the initial geometry. With increasing time, up to about 60 fm/c, the value of the v\ slope 
increases or decreases, depending on the balance between strong attractive and repulsive 
effects on particles during the rescattering process. At 40 MeV/nucleon, the net contribution 
is attractive, whereas at 150 MeV/nucleon it is repulsive, leading to the corresponding 
negative and positive flows. After 60 fm/c, the final-state interactions (FSI) still affect 
the collisions at the low beam energy 40 MeV/nucleon, and the contributions of two-body 
scatterings start to become stronger than those of the mean field potentials. It is found 
that, although the FSI affects the final flow results, the UrQMD-X can not describe well the 
data at 40 MeV/nucleon. The UrQMD-XIII calculations can describe reasonably well both 
results at 40 and 150 MeV/nucleon. 
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Finally, in Fig. [BJ the excitation function of the slope of directed flow at mid-rapidity 
for Z = 1 particles from semi-central Au+Au collisions is shown for the range of INDRA 
energies from 40 to 150 MeV/nucleon. The data are from Ref. [2| and compared to the 
results obtained with the same UrQMD-X and UrQMD-XIII whose time evolution has just 
been discussed. It is seen clearly that the UrQMD-XIII calculations describe the data well 
over the whole energy range including the location of the transition energy which is closely 
reproduced. It has to be noticed, however, that, in the energy range 90-150 MeV/nucleon 
covered also by FOPI experiments, the FOPI result for the v\ slope of Z = 1 particles is 
higher by up to 0.05 than the INDRA data j^j. Thus, in order to pin down the exact form 
of the medium modifications on both the mean-field potentials and the collision terms, it 
will also be necessary to further reduce remaining uncertainties on the experimental side 



V. CONCLUSIONS AND OUTLOOK 

In summary, by using the microscopic transport model UrQMD, we have performed a 
systematic investigation on the effects of the momentum dependence of the mean field and 
the corrections from various aspects in the collision term, especially the modification of the 
density and momentum dependence of nucleon-nucleon elastic cross sections (NNECS) on 
the observables, such as the directed and elliptic flow and the stopping in HICs at low SIS 
energies. It is clearly seen that the collective flows and the nuclear stopping are sensitive to 
all these effects, meaning that they are important to the non-equilibrium dynamic process 
in HICs at low SIS energies. 

Further, in order to describe experimental observables systematically, a consistent consid- 
eration of the uncertainties associated with both, the mean-field and the two-body collisions, 
is extremely important and should be paid more attention. The momentum dependence 
of the mean field potentials and of the density-dependent NNECS is found to be sensi- 
tive to the collectivity exhibited by the collision dynamics of HICs. At INDRA energies 
(40-150 MeV/nucleon), the dynamic transport with a soft equation-of-state with momen- 
tum dependence (SM-EoS) and with momentum- dependent and density-modified NNECS 
(cr* z (p, Pnn, 8)) describes the slope of the directed flow of Z = 1 particles at mid-rapidity 
rather well. 

From Ref. Q we learned that the cluster charge-Z dependence of the transition energy 
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of flows has been discovered. However, the verification by the theoretical side is still absent 
since in many previous theoretical investigations the vanishing flow was studied with a 
cluster charge-Z independent global quantity — the so-called directed transverse momentum 
< pf r >. In addition, the isospin effect should be visible in (the time evolution and the 
excitation function of) collective flows at low SIS energies. The knowledge of the effect of 
isospin asymmetry on both the mean field and the NNECS might be renewed under the 
new version of the UrQMD. These studies are currently underway and will be addressed in 
a forthcoming paper. 
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